Two-dimensional random walk in a bounded domain 
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In a recent Letter Ciftci and Cakmak [Eur. Phys. Lett. 87, 60003 (2009)] showed that the two 
dimensional random walk in a bounded domain, where walkers which cross the boundary return 
to a base curve near origin with deterministic rules, can produce regular patterns. Our numerical 
calculations suggest that the cumulative probability distribution function of the returning walkers 
along the base curve is a Devil's staircase, which can be explained from the mapping of these walks 
to a non-linear stochastic map. The non-trivial probability distribution function(PDF) is a universal 
feature of CCRW characterized by the fractal dimension d = 1.75(0) of the PDF bounding curve. 

PACS numbers: 05.40.Fb,02.50.-r 



Diffusion is a basic physical process [l| that moves sub- 
stances randomly from the high concentration regimes 
to the low ones. A simple mathematical model of this 
random Brownian motion 0, Q is well described by a 
random walk [3| where the walker takes unit steps suc- 
cessively in an arbitrary direction. Several variations of 
the random walk with different boundary conditions [B[ 
has been studied in o?-dimensions. Such walks can also 
be defined on a lattice where, starting from origin, the 
walker moves randomly to one of its neighbour. The indi- 
vidual random walks are known to become scale invariant 
after a large number of steps and their radial distances 
from origin follow a normal distribution. Being simple, 
the models of random walk and some of their variations, 
have found applications in several branches of science. 
Reaction-diffusion systems,percolation, network dynam- 
ics and stock fluctuations |6[ are few to name. 

Different boundary conditions [5] are known to have 
strikingly unusual effects on the longtime behaviour of 
the random walks. Recently Ciftci and Cakmak [7| 
studied the two dimensional random walk (2DRW) in 
a bounded region, where the walker deterministically re- 
turns to a pre-defined base curve near origin after cross- 
ing the domain boundary. For example, in one partic- 
ular case, the walker chooses a new coordinate (x — > 
xj y/x 2 + y 2 , y — » y/ \J 1 x 1 + y 2 ) after crossing the bound- 
ary (\x\ < 6 or \y\ < 2) and starts a fresh walk. Interest- 
ingly these walks, hereafter named as Ciftci and Cakmak 
random walk (CCRW) produce regular patterns. 

In this Letter we show that the patterns are simple 
repetition of the base curves with centers placed at all 
the points within the domain where sum of these 

integers (i + j) is even. The cumulative distribution of 
returning walkers along the base curve is found to be 
a Devil's staircase (DS) @, which could be explained 
from the mapping of CCRW to a stochastic non-linear 
map. These DS structures, which are generic features of 
CCRW, can be characterized by the fractal dimension d 
of curve which bounds the corresponding PDF. Numeri- 



cal calculations suggest that d = 1.75(0) is universal. 
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FIG. 1: Patterns : (a) uncorrelated random walk starting 
from the base curve (b) corresponding CCRW after 10 
steps. 

In CCRW, the walker initiate a discrete time walk on 
a two dimensional square lattice starting from origin, by 
taking a random unit step ±1 both in x and y direction. 
Thus in each time step its coordinate changes from (x, y) 
to (x + a x ,y + cr y ), where a XiV = ±1 chosen randomly. 
The walk is confined in a bounded domain, usually taken 
as rectangle \x\ < b x and \y\ < b y . If the walker crosses 
the boundary , e.g. when \x\ > b x or \y\ > b y , it returns 
immediately to a new coordinate (f(6),g(9)), where 9 = 
tan~ 1 (y/x). The new co-ordinates define a curve 

x = f(6); y = g(9) (1) 

parametrized by 9, which will be referred to as the base 
curve (BC). Case-I of [7| corresponds to f(9) = sm9 
and g(9) = cos9 which generate regular patterns in a 
bounded domain, compared to simple two dimensional 
random walk which visit only a set of integer points {v} = 
(x — i,y = j) where (i + j) is even. 

The walker in CCRW eventually returns to the base 
curve and starts a fresh walk from there. Thus, it is natu- 
ral to expect that the patterns of CCRW can also be gen- 
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erated by a set of uncorrelated random walk (URW) of a 
large number of walkers starting from the base curve. In 
URW, however, one must use a absorbing boundary con- 
dition that the walk terminates when the walker crosses 
the boundary. Starting from any arbitrary point r on 
the base curve each walker visits only the points {f+ v} 
Hence, URW produces a pattern where the base curve is 
shifted by the vectors {€}. This pattern is independent 
of the initial distribution of walkers on the base curve. 

To show that URW can generate any desired pattern 
we take an example, 



Boundary 
Base curve 



|x| < 3 and \y\ < 2 
f{9) = asin(0)cos 2 (0) 
g{6) = asin 2 (0) cos(0). 



(2) 



Let the distribution of walkers P{9) along the base 
curve at an angle is be uniform [i.e 9 is a random num- 
ber uniform in (0, 2n)]. The patterns generated by this 
URW with a = 1 is shown in Fig. [TJ In this pattern the 
repeating base curves are non-overlapping as their lobes 
are bounded within a circle of radius l/v2. One can, 
in fact, generate more complicated and overlapping pat- 
terns either by choosing a > y2 or by taking different 
BCs. 
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FIG. 2: Distributions of returning walkers P{6) = P(—0) 
along the BC ([2]), and a unit circle are shown in (a) and (b) 
respectively. Clearly, P(0) is symmetric about 9 — tv/2. The 
arrows in (a) point to some forbidden regions where P{9) — 0. 
Here, PDFs are calculated with Ad = 2tt10~ 5 and N = 10 s . 
The insets of (a) and (b) show variation of (j> with N for 
respective BCs (e = 10 -4 ). 

Patterns similar to Fig. 1(a) can also be produced by 
CCRW. The walker in this case initiates the walk from 
origin and immediately after crossing the boundary it re- 
turns to the base curve The return is deterministic 
as 9 = tan _1 (y/x) in Eq. @ depends on the final co- 
ordinates (x, y) of the walker. Corresponding pattern is 
shown in Fig 1(b). Unlike Fig 1(a), here, the patterns ap- 
pear discontinuous, which indicates that certain regions 
along the BC are never visited. This is also confirmed 
from the numerical results of the distribution of return- 
ing walkers P(9) along the BC in the range 9 =€ (0, tt). 
Fig. [5] shows the distribution P(9), where some of the 
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FIG. 3: P(8\A0) does not approach any definite limit as 
A9 — s> 0. (a) and (b)correspond to CCRW with base curve 
p| and a unit circle respectively. 



forbidden regions [P(9) = 0] are marked with an arrows. 
A quantitative measure of the forbidden regions is 
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Q(e-P(6))d0, 



(3) 



where e ~ is a pre-determined positive number and 
0(x) is the Heaviside step function. In fact <j> saturates 
to cj) s = 0.232(4) as the number of time steps N increases 
[Inset of Fig. Ufa)]. A non-zero <j> 3 is not a typical char- 
acteristic of CCRW. When the BC is a circle (Case -I 
of 0] ) , we find that the walker returns to all values of 9 
with finite probability [as shown in Fig. 2(b)] and thus 
<t> s -> [Inset of Fig. 2(b)]. 

A unique feature of CCRW that emerges from the dis- 
tribution of returning walkers shown in Fig. [5] is that 
this probability measure may not be represented in any 
functional form. To get a functional form, namely the 
PDF, one needs to count the fraction of walkers P(0|A0) 
that comes back to a bin of size A9 about 9 and then 
take the limit [H| 



P{9) 



lim P(0|A0). 

A0^O 



(4) 



The PDF is well defined only when above limit exists. In 
Fig. 3(a) we have plotted P(0| A0) against A0 for CCRW 
with BC © and for two different values of 9 = 0.7, 0.4. 
It is evident that the limit A0 — > does not exist. Figure 
3(b) there shows the same for Case-I of Q and for 

0.7. ().:'>. 

It is natural to ask, if the cumulative distribution de- 
fined by P(9 > (j>) = J" P{9)d9 is a well defined function. 
Fig. 4 shows that P{9 > 4>) against <j> for the base curve 
/(0) = sin{9),g{9) = cos(9) is a Devils staircase [§], 
not differentiable at infinitely many points. Thus it is 
not surprising that P(0|A0) does not have a well defined 
limit for A0 0. 

Contrary to naive expectation that any deterministic 
dynamics, when added to an existing stochastic system, 
does not alter the the stochastic behavior, here we ob- 
serve that the deterministic return of a two dimensional 
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FIG. 4: Cumulative distribution P(9 > (/>) for base curve 
f(8) = cos(8) and g(6) = sin(#) is a Devil's staircase. Insets 
show enlarged part of different regions. 



random walker to a pre-defined base curve partially de- 
stroy the random-ness resulting in a non-trivial probabil- 
ity distribution. To understand it better, let us take the 
following CCRW, 



Boundary 
Base curve 



\x 

X 



< 1/2 and \y\ < 1/2 
= cos(9) /2 and y = sin(6)/2 (5) 



where 9 is angle made with x-axis by the returning ran- 
dom walker. In this case, the walker crosses the bound- 
ary in every attempted walk and returns to the base 
curve. Effectively the walker traces different points on 
the base curve using a stochastic map Xt+\ — (xt + 
°f )l r u Vt+i = (yt + Vt)l r t where a x,v = ± chosen ran- 
domly (by the two dimensional walker during its walk) 
and r t — \J (x t + erf) 2 + (y t + of). One can, in fact 
write an equivalent one dimensional map in terms of 
2 = y/x = tan(0) as 



with h(z, a x , a v ) 



Z t+ 1 = fc(z t , of ,pf ) 

z + 2crVl + z 2 



1 + 2a x Vl + z 2 



(6) 



On each iteration z takes a new value by choosing one of 
the four non-linear functions h{z, ±1, ±1) randomly as 
shown in Fig. [SJ These functions h(z,<r x ,a v ) have at- 
tractive fixed points at z* = a x a y . Note, that if a x were 
not random (say a x = 1) then z would evolve using the 
map z t +\ — h(zt,^,Ot) and ultimately remains confined 
in the region \z\ < 1 as t — > oo. Addition of the other 
two functions h(z, —1, a v ) allows z to move out of the re- 
gion \z\ < 1. Corresponding distribution of 6 = tan^ 1 {z) 
(shown as a cumulative distribution in Fig. [5]) compares 
well with Fig. HI which indicates that CCRW © has all 
the characteristic features of the generic two dimensional 
bounded random walk. 

We argue that this non-trivial behavior, which occurs 
even for CCRW ([5]) , is the effect of the non-linearity that 
exists in the stochastic map (j6]). In the following we 
show that a linear stochastic map having identical fixed 
points would not show any such features. In particular, 
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FIG. 5: The stochastic non-linear map Q corresponding to 
CCRW ©. The fixed points of h(z,a x ,a y ) are a x a y . 



the steady state probability of z has a well-defined PDF. 
Let us take the following stochastic linear map 



(7) 



which has the same fixed points as (J5]) and calculate the 
steady state distribution g(z). Direct iteration of this 
map yields 



z t +i 



t 

E 

71 = 



<j x 1 

2 '2 t - n 



n 



fe=n+l 

Since product of a x s produces a random sign ±, we have 



Zt+l 



/ < On v 



n=0 



T 

n 
2 



(8) 



where the term containing the initial value Zq (being ex- 
ponentially small for large t) is dropped, and we have 
taken r x,y = <J x ^ n . The right hand side of the above 
equation resembles the Hamiltonian of a spin system on 
a lattice of size L = t, where the spins t x ' v of two differ- 
ent kinds of particles x and y at the site n interact with 
an inhomogeneous magnetic field B n = 2~ n . Thus the 
expression for the energy of the system is 



£ = 5>rf(l-^). 

i=0 



(9) 



The energy distribution of this model P(E; j3) in equi- 
librium is related to the stationary distribution g(z) of 



j(z) = lim P{E = z;P). 

,8-KI 



(10) 



Since, for classical systems with partition function Zl{P), 



P{E,P)=C~ 1 



Z L (f3) 



E 



where C 1 [f(s);x] = f°°e sx f(s)ds is the inverse 
Laplace transform, we have 



g(z) = c- 



Zl{s) 
Zl(0) 



(ii) 
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FIG. 6: Cumulative distribution in the steady sate of Eq. ([6]) 
which is equivalent to CCRW (|5|, is compared with Eq. 



The partition function of the model (j9)) 

sinh/3sinh2/3 



Z L {fi) 



sinh(2- L /3) sinh(2- L 2^) 
is used in Eq. (fTl"]) to obtain 



9(z) = 



1/4 

) 3-e(z)z 
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|z|<l 
1< \z\ < 3 

kl >3, 



(12) 



in thermodynamic limit L oo. Here e(z) is the 
signum function. Correspondingly, the distribution of 
6 = tan -1 z is P(&) — g(tan#) sec 2 9, and its cumulative 
distribution is 



2+tanQ) 

P{6 > (j)) = I 8+6 tan 0+e(tan <f>) sec <j> 2 





Hi 



|tan0| < 1 

1 < | tan (f>\ < 3 

I tan 6\ > 3 



In Fig. [6] we have plotted the cumulative distribu- 
tion P(# > 4>) for both the non-linear map (j6]) and the 
stochastic linear map (|8|) which has identical fixed points. 
It is quite evident from the figure that P(6 > (f>) for the 
non-linear map, which corresponds to the CCRW ([5]), has 
a structure of Devil's staircase, whereas the same for the 
linear map is a continuous and differentiable function. 
Thus, it is suggestive that the non-trivial distribution of 
returning walkers on the base curve is an artifact of non- 
linearity. 

To characterize the nontrivial distribution P(6) of 
CCRW we use the fractal dimension of its bounding 
curve. From naive box-counting methods one expects 
that the bounding curve, which is a fractal, is covered by 
v ~ A<fi~ d segments of size A<fi with d / 1. For CCRW 
with BC a plot of v versus A(j> in log-scale (Fig. 6) 
shows that d = 1.75(8). Same calculation for a few other 
BCs [|[ results in d ~ 1.75, which made us to conjec- 
ture that possibly the fractal dimension of PDF bounding 
curve of CCRW is universal (d = 7/4). 



In conclusion, we explain that patterns of Ciftci- 
Camak random walk [7] (CCRW) are simple repetition of 




o.ooi 



A(j)/27t 



FIG. 7: Plot of v versus A0 in log scale for the PDF of CCRW 
with BC ([2]), a unit circle, and a unit square Q are shown 
with symbols x , o and □ respectively. A line with slope —7/4 
is drawn to guide the eye. 



the base curve. An unique feature of CCRWs is that the 
distribution of returning walkers along the base curve can 
not be represented by a functional form; corresponding 
cumulative distributions are Devil's staircases. A cor- 
respondence of CCRW with stochastic Wiley non-linear 
maps reveals that this unusual distribution is an artifact 
of non-linearity. A quantitative characterization of the 
distribution could be fractal dimension d of its bound- 
ing curve. Our numerical calculations of CCRWs with 
different base curves show that d ~ 1.75. 
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